Spin-distribution functionals and correlation energy of the Heisenberg model 
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We analyse the ground-state energy and correlation energy of the Heisenberg model as a function 
of spin, both in the ferromagnetic and in the antiferromagnetic case, and in one, two and three 
dimensions. First, we present a comparative analysis of known expressions for the ground-state 
energy Eo(S) of homogeneous Heisenberg models. In the one-dimensional antiferromagnetic case 
we propose an improved expression for Eq(S), which takes into account Bethe-Ansatz data for 
S = 1/2. Next, we consider inhomogeneous Heisenberg models (e.g., exposed to spatially varying 
external fields). We prove a Hohenberg-Kohn-like theorem stating that in this case the ground-state 
energy is a functional of the spin distribution, and that this distribution encapsulates the entire 
physics of the system, regardless of the external fields. Building on this theorem, we then pro- 
pose a local-density-type approximation that allows to utilize the results obtained for homogeneous 
systems also in inhomogeneous situations. We conjecture a scaling law for the dependence of the 
correlation functional on dimensionality, which is well satisfied by existing numerical data. Finally, 
we investigate the importance of the spin-correlation energy by comparing results obtained with the 
proposed correlation functional to ones from an uncorrelated mean-field calculation, taking as our 
example a linear spin-density wave state. 

PACS numbers: 71.15.Mb, 75.10.Jm, 75.50.Ee, 75.30. Fv 



I. INTRODUCTION 

In this paper we study the ground-state energy, cor- 
relation energy and related quantities of the Heisenberg 
model. The homogeneous Heisenberg model is defined by 
the Hamiltonian 



The inhomogeneous Heisenberg model, characterized 
by broken translational symmetry, is obtained by adding 
a spatially varying magnetic field to Hq, 



H 



(2) 



H 



(ij) 



(1) 



where the Si are spin vector operators satisfying 
Sf\Sm) — S(S + l)\Sm) and Si <z \Sm) — m\Sm), and 
S and m are the spin quantum numbers of the parti- 
cles under study. Although in accordance with common 
terminology the Si are called spin operators, they re- 
ally represent total angular momentum and are not re- 
stricted to be of purely spin origin, (ij) indicates a sum 
over nearest neighbors on a lattice of dimensionality d, 
and J is the spin-spin interaction constant, parametrizing 
the exchange interaction of the underlying microscopic 
Hamiltonian. 1-3 For antiferromagnetism J > 0, while 
for ferromagnetism J < 0. This model was originally 
proposed in 1926 to explain ferromagnetism in transi- 
tion metals, 4,5 but has since then found a large num- 
ber of other applications to the magnetic properties of 
matter. 1-3 Recent examples are antiferromagnetic chains 
in complex oxides and other low-dimensional magnets 6 
or studies of magnetic effects on crystal-field splittings in 
rare-earth compounds. 7 



where Bi can either be an externally applied field or an 
internal field due to magnetism in the system. A vari- 
ety of different sources and manifestations of the inhomo- 
geneity Bi has been studied in the recent literature, often 
in conjunction with the synthesis and investigation of real 
materials, whose magnetic properties are necessarily spa- 
tially inhomogeneous, but still to some extent describable 
by modified Heisenberg models. 6,8 The interpretation of 
experimental data for realistic materials in terms of the 
Heisenberg model, in particular in the presence of stag- 
gered or otherwise spatially varying magnetic fields, must 
be based on a solid understanding of the behaviour of the 
model (2) in the presence of inhomogeneity. 

The homogeneous Heisenberg model in d = 1 dimen- 
sion and for S = 1/2 has an exact analytical solution in 
terms of the Bethe Ansatz, 9,10 but the same Ansatz does 
not work in higher dimensions, in which no exact solution 
is known. It is also hard to generalize to inhomogeneous 
situations. In the present paper we combine exact and 
approximate results obtained within a variety of differ- 
ent approaches and techniques, to provide a systematic 
analysis of the ground-state and correlation energy of the 
Heisenberg model both in the ferro and in the antiferro- 
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magnetic case and for d = 1,2,3 dimensions. 

In Sec. II we provide a comparative analysis of avail- 
able expressions for the ground-state energy of the homo- 
geneous Heisenbcrg model as a function of spin S, dime- 
nionality d, and coupling constant J. This section also 
contains a proposal for an improved expression that goes 
beyond those available in the literature in taking into ac- 
count a Bethe Ansatz result for S = 1/2. In Sec. Ill we 
use concepts of density-functional theory to extend the 
utility of the homogeneous results discussed in the pro- 
ceeding section to inhomogcneous situations. This sec- 
tion includes a proof of a Hohcnbcrg-Kohn-type theorem 
for a large class of generalized Heisenberg models, and 
a proposal for a simple local-density approximation. In 
Sec. IV we then study the correlation energy in homoge- 
neous and inhomogeneous Heisenberg models, in order to 
assess the importance of correlations and quantum fluc- 
tuations as a function of dimensionality and spin. As 
an explicit example for a physically interesting type of 
inhomogcneity we consider a linear spin-density wave, 
and explore the differences between a mean-field and a 
density-functional treatment of the resulting inhomoge- 
neous spin distribution. 

II. HOMOGENEOUS HEISENBERG MODELS: 
GROUND-STATE ENERGY AS FUNCTION OF 
SPIN 

This section provides a brief review of what is known 
about the ground-state energy of homogeneous Heisen- 
berg models. Although the expressions collected below 
for J > and J < and d = 1,2,3 are given at various 
places in the literature, we have not found a systematic 
collection and comparison at one single place. For the 
convenience of the reader and future reference we there- 
fore provide such a comparison below. We also add a 
new expression to the list [Eq. (15)], which is a slight 
improvement on one of the earlier results. 

A. Ferromagnetic case 

Let us first consider the ferromagnetic case, which is 
much simpler than the antiferromagnetic one. At zero 
temperature all spins are parallel, and the correspond- 
ing spin operators commute with each other. The exact 
ground-state energy is then the same as that obtained in 
the mean-field approximation, and is given by 

E™{S,J,z) = J Z -NS\ (3) 

where z is the number of nearest neighbours, which on 
linear, square and cubic lattices is related to the dimen- 
sionality by d = z/2 and, as above, J < for ferromag- 
netism. In Fig. 1 we show the resulting curves for the 
energy per site and interaction strength 



4 M {S, d) = J-E FM (S, J, z = 2d) = -dS\ (4) 
for one, two and three dimensions. 

B. Antiferromagnetic case 

The antiferromagnetic (AFM) case is much more com- 
plicated than the ferromagnetic one, because in spite 
of its name the ground-state of the antiferromagnetic 
Heisenberg model is not simply the 'antiferromagnetic' 
state consisting of alternating spin up and spin down 
states with respect to a fixed direction (i.e., the Neel 
state), but a quantum superposition of states involving 
also spins along the perpendicular axes. In the one- 
dimensional S = 1/2 case the structure of the corre- 
sponding ground-state for N — > oo is known exactly 
by means of the Bethe Ansatz. 3,9 ' 10 The corresponding 
ground-state energy is 

E AFM (S= 1 -,J,d=l) 

= E FM (S = ^, J, z = 2) - JN In 2. (5) 

For the energy per site and interaction strength one ob- 
tains from this 10 

er M (S = \,d=l) = ^jE^iS =\,J,d=l) (6) 

= - - In 2 = -0.44314718. (7) 

4 

Similar exact results are not available in the general 
case, for arbitrary S and d. However, a set of useful ap- 
proximate expressions was derived in two early papers by 
Anderson. 11 ' 12 In the first of these it is shown by means 
of a variational argument that the energy of the antifer- 
romagnetic ground state must lie in the interval 11 

-^(l + ^<er M (S,z)<-^S\ (8) 

where z — 2d is again the coordination number of the 
linear, square or cubic lattice under study. A simple es- 
timate is obtained by using the center of this interval, 12 
i.e., 

but the quality of this estimate, which is quite good for 
d = 3, deteriorates for d = 2 and d = 1. A numerical 
calculation based on spin-wave theory leads to the more 
precise results 12 

e$ FM (S, d=l) = -S 2 +(--l)s (10) 
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= -5 2 1 + 



0.36338 



e£ FM (S,d = 2) 
e£ FM (S,d = 3) 



-2S 2 ( 1 + 
-35 2 ( 1 + 



5 
0.158 



5* 
0.097 
5 



(11) 
(12) 

(13) 



We can assess the quality of these expressions by sub- 
stituting 5 = 1/2 in the first of them. The result 



S AFM 
-0 



(S 



Id 



1 



-0.43169 is within 2.6% of the 



exact Bethe Ansatz value reported above. 

Indeed, significant improvement over Eq. (10) has only 
been obtained recently, with the use of modern comput- 
ing facilities and advanced numerical techniques. In such 
work Lou ct al. 13 used (50 years after Ref. 12) the density- 
matrix renormalization group (DMRG) to calculate cor- 
rections to Eq. (10) for values of 5 ranging from 1/2 to 
5 in steps of 1/2. These authors propose a fit to their 
numerical data which in our present notation reads 



So FM ' fitl (S, d = 1) = -S' J + ( ^ - 1 ) 5 - 0.03262 



-0.0030-j; - [0.338 - ^ ~ 



e~ TrS cos(2 7 rS*). (14) 



We note that for 5 



AFMJitl 

-o 



(5 



,d = 1) 



1/2, where it predicts 
-0.516459, this fit is actu- 



ally worse than the earlier expression (11), deviating by 
17% from the exact Bethe Ansatz value. On the other 
hand, as shown in Ref. 13 by comparison with DMRG 
data and other highly precise numerical results, the fit is 
excellent for higher values of 5. 

In order to obtain a closed expression that can also be 
applied at 5 = 1/2, a slight modification to the fit by Lou 
et al. is sufficient. To this end we propose the alternative 
expression 



AFMJitl 
3 



(S,d = l) = -S 2 + 



15 



-0.03262 - 0.0030-^ 



0.0015-^3 



„ 000 0.28 0.035 
0.338 - — + -^ 



e -7rS cos(27rS), 



(15) 



which differs from (14) in the inclusion of two cubic terms 
in 1/5. The value at 5 = 1/2 predicted by this expres- 



AFMJU2 

3 



(S 



,d = 1) = -0.446253, deviates by 



only 0.7% from the exact Bethe Ansatz value. 

Figures 2 and 3 display the various AFM energy ex- 
pressions collected above. In Fig. 2 we compare the rigor- 
ous but rather wide interval provided by expression (8), 
with the spin-wave results (11) - (13) and (in d = 1) 
the DMRG data of Ref. 13. Obviously, on this scale the 
spin- wave expression (11) already provides an excellent 
approximation to the DMRG data. Neither of the two fits 
(14) and (15) can be distinguished from the much simpler 
expression (11) on the scale of the figure. Unfortunately, 



no highly precise numerical reference data, similar to the 
DMRG results of Ref. 13, seem to be available in d = 2 
and d — 3, but the approximations leading to the sim- 
ple analytical formulae given above are expected to work 
better as d increases. This expectation is corroborated 
by noting that the rigorous interval (8) shrinks with in- 
creasing d. 

Interestingly, both the numerically highly precise 
DMRG data of Ref. 13 (in d = 1) and the spin-wave 
expressions (11) to (13) (in d = 1, 2, 3) systematically lie 
closer to the more negative boundary of the interval than 
to the less negative one, and in d = 1 the DMRG values 
are still a little closer to this boundary than the curve 
predicted by Eq. (11). This shows that the lower bound 
in Eq. (8) is tighter than the upper one. The simple esti- 
mate (9), on the other hand, by construction falls in the 
middle of the interval and becomes less reliable for lower 
d. In the interest of readability we have not displayed 
the curves corresponding to this estimate in the figures. 

Figure 3 shows that on a smaller scale the differences 
between the more precise expressions become important. 
Here we compare the two fits, (14) and (15), to the 
DMRG data. To make the details of the fits, and the 
interesting oscillatory structure they display, clearly vis- 
ible, we have subtracted the spin-wave expression (11), 
which is common to both fits. The fit proposed above, 
Eq. (15), is slightly inferior to the one developed by Lou 
et al., Eq. (14), around 5 = 3/2, but unlike the latter 
recovers the exactly known 5=1/2 data point to within 
less than one percent. 

In the present paper we are mainly concerned with the 
homogeneous or inhomogencous Hcisenberg model on lin- 
ear, square and cubic lattices. Expressions (or numerical 
values) for the ground-state energy can also be derived 
for many other variations of the Heisenberg model, such 
as lattices with helical boundary conditions, 14 or with 
anisotropic interactions. 15 Although we do not consider 
such models in the present paper, many of our results 
can be extended to them in a straightforward way. 



III. INHOMOGENEOUS HEISENBERG MODELS: 
SPIN-DISTRIBUTION FUNCTIONALS 

Based on the analysis of the preceeding section we rec- 
ommend the use of expressions (11) to (13) in calcula- 
tions requiring simple expressions for the ground-state 
energy of the homogeneous antiferromagnetic Heisenberg 
model in one, two, and three dimensions. In one dimen- 
sion, where the simple expressions fare worst, either of 
the two fits (14) and (15) provides a significant improve- 
ment in accuracy, but only the latter recovers the Bethe 
Ansatz value at 5 = 1/2. However, the utility of any 
of these expressions is rather limited due to the restric- 
tion to spatial homogeneity. Externally applied magnetic 
fields that vary in space, calculations of magnetic effects 
on crystal-field splitting, description of nontrivial inter- 
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nal order, etc., require use of the inhomogeneous Heisen- 
berg model. 8 Unfortunately, if translational invariance is 
broken the Bethe Ansatz, spin-wave theory, DMRG and 
most other approaches encounter very significant compu- 
tational difficulties. 

In the case of ab initio calculations a many-body tech- 
nique that has had considerable success in the application 
to inhomogeneous systems is density-functional theory 
(DFT), 16-18 but it is not very common to apply DFT 
also to model Hamiltonians. However, following pioneer- 
ing work by Gunnarsson and Schonhammer, 19 DFT was 
recently formulated and applied for the one-dimensional 
Hubbard model. 20-22 In this section we build on this ex- 
perience to explore how DFT can become a useful tool 
also in studies of the inhomogeneous Heisenberg model. 
To this end we prove, in subsection III A, a Hohenberg- 
Kohn-type theorem for a wide class of generalized Heisen- 
berg models (of which the models discussed above are 
special cases). In subsection IIIB we then use this the- 
orem and the explicit expressions discussed in Sec. II to 
construct a simple local-density approximation for inho- 
mogeneous Heisenberg models. 



A. Hohenberg-Kohn theorem for generalized 
Heisenberg models 

A first question that must be answered before DFT can 
be usefully employed is what the fundamental variable is. 
In ab initio DFT one mostly chooses the particle density 
n(r) or its spin- resolved counterpart n CT (r), 16-18 although 
other choices are occasionally useful. 23-25 In the case of 
the Hubbard model the basic variable is the site occu- 
pation number rij. 19-22 In the present case we propose 
to use the spin vector Si, which is the only fundamen- 
tal dynamical variable appearing in the definition of the 
Heisenberg model. 

In the interest of generality, in the present section we 
consider a generalized Heisenberg model of the form 



H g — ^2 Jij&i • Sj + Si 



B, 



(16) 



Unlike in Eqs. (1) and (2) the sum in the first term on the 
right-hand side is not restricted to nearest neighbours, 
and the interaction can depend in any way on the 
indices of the involved sites. In particular, it can extend 
to next-nearest neighbours and beyond, or alternate be- 
tween ferromagnetic and antiferromagnetic along some 
direction in the crystal. Both of these features arc found 
in realistic magnetic crystals. This relaxation of con- 
straints on J may appear a considerable complication, 
but it turns out that the proof of the Hohenberg-Kohn 
theorem is essentially unaffected by the extra generality. 
(As pointed out above, we do not consider anisotropic 
Heisenberg models, in which J couples differently to dif- 
ferent components of S, in this paper, but the general- 
ization of the theorem to this case is straightforward.) 



Following the steps of Hohcnberg and Kohn, we now 
consider two Hamiltonians with same interaction , but 
exposed to two different magnetic fields Bi and BJ. Thus 



^2 J ijSi ■ S J ■ + ^2 Si ' Bi 



* = £ * (17) 



*' = E' *', (18) 



where Eq and E' are the ground-state energies in the 
fields Bi and BJ, and * and are the corresponding 
ground-state wave functions. 

As a consequence of the variational principle we have 
the inequality 



E = < <*'|ff|tf'}, 



(19) 



since *' is not the ground-state wave function belong- 
ing to H (assumed nondegenerate) . By adding and sub- 
tracting the term J2i S» ■ B • on the right-hand side, the 
inequality becomes 

E = (*|#|*) < <tf'|iJ'|tf') + 2<*'|Si • (Bj - B{)|*'>. 

i 

(20) 

Here the first term on the right-hand side is just the 
ground-state energy Eq, of Hamiltonian H'. With the 
abbreviations AB; = B, - BJ, S» = (*|Si|*) and SJ = 
("J' | Si I*'} the preceeding equation then becomes 



E < E' + J2 ^ ■ AB,. 



(21) 



Now we repeat the same argument starting with the 
Hamiltonian H' . The variational principle guarantees 
that 



E' Q = (ty'\H'\V') < (*|^'|*). 



(22) 



By adding and subtracting the term J2 i Sj-Bj, we obtain, 
in the same way as before, 

E' a < E -J2 S * ■ AB - ( 23 ) 

i 

Addition of Eqs. (21) and (23) leads to 

E + E <S + B + ^(SJ-S J )AB J . (24) 



If we now assume that SJ = Si, i.e., that the two spin dis- 
tributions corresponding to the two different wave func- 
tions ^ and *' are identical, then the previous equation 
reduces to the contradiction 



Eq + Eq < Eq + E' Q . 



(25) 
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This contradiction shows that two distinct nondegen- 
erate ground states can never lead to the same spin distri- 
bution. Hence, given some arbitrary spin distribution Sj 
there is at most one nondegenerate wave function which 
gives rise to it. In other words: the spin distribution 
uniquely determines the wave function. This means that 
the wave function is a functional 26 of the spin distri- 
bution, i.e., "J — "J [S^. This is the statement of the 
Hohcnberg-Kohn theorem for the Heisenberg model. 

For completeness we mention that the above proof by 
contradiction, patterned after the one first presented by 
Hohcnbcrg and Kohn in the ab initio case, 27 is not the 
only possible one. The constrained-search technique of 
Levy 28 and Lieb 29 is also easily adapted to the present 
case. The ground-state wave function in this approach is 
uniquely defined by its spin distribution as the wave func- 
tion that minimizes (^\H\^) and reproduces Sj. This 
minimization defines the functional 

E LL [Si]=mm{*\H\*). (26) 

whose minimum is the ground-state energy. 

An immediate consequence of either formulation of the 
proof is that the ground-state expectation value of any 
observable O, is also a functional of the spin distribution, 
defined via 

0[St] = (*[S;]|6|*[S;]), (27) 

and this functional is the same regardless of the strength 
and direction of the magnetic field Bi, i.e., it is universal 
with respect to external fields. Note that the theorem 
applies to any ground-state observable. For example, it 
implies that also all multi-spin correlation functions of 
the general form 

Cn,n+l,n+2,... '■= | S n S„+i S„+2 • • • |^) (28) 

are uniquely determined by the single-spin expectation 
value S„ = (*|S n |*). This is trivially true for the 
nearest-neighbor correlation function C^°^_ 1 of a homo- 
geneous one-dimensional system, which as a consequence 
of the definition of the homogeneous Hamiltonian (1) is 
simply given by 

CtT+i - j^Eo(S), (29) 

but the above proof guarantees that the more compli- 
cated correlation functions involving more than two spins 
and/or spatially inhomogeneous spin distributions are, in 
principle, also functions of Sj only. 

B. Local-density approximation 

Another consequence of this Hciscnberg-model for- 
mulation of the Hohenberg-Kohn theorem is that the 
model's ground-state energy and spin distribution can 



be obtained by application of the variational principle 
to spin distributions, instead of wave functions. In com- 
plex situations this can be a major simplification, but 
to extract this information is, of course, still highly non- 
trivial. The most straightforward thing to do would be 
to set up an approximation for the total energy of the 
system under study as a functional of the spin distribu- 
tion, and to minimize with respect to Sj. In ab initio 
DFT this is not the preferred way to proceed because it 
turns out to be hard to conceive good density function- 
al for the kinetic energy. In practical applications of ab 
initio DFT one therefore commonly employs an indirect 
minimization scheme, leading to the widely used Kohn- 
Sham equations. Although this could also be done in the 
present case, there is no need for introducing a Kohn- 
Sham system for the Heisenberg model, since there is no 
kinetic-energy term in the first place. Direct minimiza- 
tion of total-energy expressions seems much more conve- 
nient (and more analogous to the way the model is usu- 
ally treated in statistical physics) than indirect minimiza- 
tions. In order to explore how such a direct minimization 
can proceed we first construct, in this subsection, a local- 
density approximation for the simpler Heisenberg models 
discussed in the preceeding section. 

Let us write the two contributions to the ground-state 
energy of the inhomogeneous Heisenberg model as Ej 
and Eb, which arc the ground-state expectation values 
of the first and second term on the right-hand side of 
Eq. (2), respectively. The mean-field approximation for 
Ej yields 

E** F [Si] = jJ^Si-Sj + Y^Si-Vi 

ij i 

= Ey F [Si] + E B [Si], (30) 

where, as above, Si = (\I/|Sj|\I'). In Sec. IV we quantify 
the error made by the neglect of correlation effects arising 
from use of Ef F [Si] in place of Ej [Si] . 

A guideline for the construction of better functionals 
than Sj fF [Si] is provided by ab initio DFT 16 or recent 
work on the Hubbard model. 20-22 In the former, the total 
energy of an arbitrarily inhomogeneous system is written 
as 

E [n{r)} = T s [n(r)} + E H [n(r)} + E v [n(r)} + E xc [n(r)}, 

(31) 

where T s is the noninteracting kinetic energy, 

BrtM-\j*rj*S?8&i (32) 

the Hartree energy, and E v the potential energy arising 
from the external field v(r). The local-density approxi- 
mation (LDA) for the exchange-correlation (xc) energy 
is 

E xc [n(r)} » E^ A [n{v)} = J d 3 r e xc (n)\ n ^ n(r) . (33) 



•5 



This expression locally substitutes the xc energy of the 
inhomogeneous system by the one of a homogeneous sys- 
tem of same density. The necessary input expression for 
e xc (n) is obtained by subtracting the Hartree and nonin- 
teracting kinetic energy from the ground-state energy of 
the homogeneous system, eo(n). 

In the present context we write, in analogy to Eq. (31), 



E, 



[S i ]=E? F [S i ]+E B [S i ]+E c [S i 



(34) 



where E ] f F is defined in Eq. (30), Eb is the potential en- 
ergy arising from the external field , and E c is by def- 
inition the difference between the mean-field result and 
the correct one, i.e., the correlation energy. (There is no 
Heisenberg-model counterpart to the kinetic energy term, 
and we avoid the expression 'exchange-correlation en- 
ergy' because in common terminology the entire Heisen- 
berg Hamiltonian is due to 'exchange'.) 

To obtain an explicit scheme we now propose to ap- 
proximate, in analogy to (33), 



i? c [S,]«^[S^^e c (5)|^| Ss |, 



(35) 



where e c (S) is obtained by subtracting the mean- field en- 
ergy, —dS 2 , from the homogeneous expressions for eo(5) 
discussed in Sec. II. As an explicit example, the LDA 
approximation for the correlation energy of an inhomo- 
geneous antifcrromagnetic Heisenberg model in one di- 
mension becomes 



E c ,AFM,d=l[Si] — J 



1 



EiSii, 



(36) 



where we used Eq. (10) for 60(5*). Of course Eqs. (14) or 
(15) can be used in the same way in d = 1, and Eqs. (12) 
and (13) in d = 2 and d = 3, respectively. The full 
ground-state energy is then for any d approximated as 

£ [S t ] « £ iDA [S t ] = E™ F [S i ] + E B [S t ] + EZ DA [Si\. 

(37) 

Clearly Eq. (35) is a rather simple approximation, 
whose quality may vary widely depending on the circum- 
stances (e.g., values of J and d, or spatial dependence 
of Bj). At present it is motivated mainly by the con- 
siderable practical success of its counterpart in ab inito 
DFT, 16-18 and by the encouraging results obtained re- 
cently with a Bethe-Ansatz based LDA for inhomoge- 
neous Hubbard models. 20-22 It is clear, however, that the 
LDA contains essential correlation effects not accounted 
for by the mean-field expression (30). In spite of the 
extra term, minimization of (34) with (35) is no more 
complicated than that of (30). Eq. (35) thus shows one 
way in which the expressions listed in Sec. II for homo- 
geneous systems can be applied to inhomogeneous situa- 
tions. A simple application of these ideas is worked out 
in Sec. IV B. 



C. A scaling hypothesis 

An interesting feature of the functional obtained by 
combining Eq. (35) with the explicit formulae of Sec. II 
is that the resulting expressions depend explicitly on the 
interaction J and the dimensionality d. The dependence 
of the ab initio functionals on these parameters is not well 
known, and in particular the ^-dependence of the func- 
tional is still subject of many ongoing investigations. 30 
Even in the much simpler case of the Hubbard model 
it is only the interaction dependence which is featured 
explicitly in the available functionals, 19-22 whereas the 
dependence on dimensionality is essentially unknown. In 
this context it may be useful to have, for the Heisen- 
berg model, a number of simple expressions that depend 
explicitly on dimensionality and interaction, so that the 
role of these parameters in the functional can be explored 
in a simplified environment. 

As an explicit example, we consider scaling properties 
of the functional as a function of dimensionality d. The 
Hartree-like term oc S 2 in Eqs. (11) to (13) clearly scales 
linearly with d. Hence 



ey F (d)=de<r(d=l) 



MF, 



(38) 



Interestingly, the correlation energy contribution cx S 
also obeys a similar, albeit less obvious, scaling law. 
From the explicit expressions (11) to (13) for eQ FM (d) 
one obtains for the relation between the ratio of the cor- 
relation energies and the ratio of the dimensionalities 



e c (d = 2) 
e c (d = 1) 



= 0.870 = 



and 



e c (d = 3) /3 
— -f = 0.921 = - 
e c (d=2) \2 



(39) 



(40) 



where x and y are exponents to be determined. Nu- 
merically one finds x = —0.201 and y = —0.203. The 
near-equality of these two exponents among each other 
and to the integer fraction —1/5 leads us to conjecture 
the following dimensional scaling law: 



e c {d) = d- r >e c {d=l), 



(41) 



where the scaling exponent n = 1/5. This scaling law 
accounts for the numbers in Eqs. (11) to (13) to within 
w lO- 3 . 

Of course, at present the scaling law (41) is only a 
conjecture, but one that is consistent with the num- 
bers of spin-wave theory. It also correctly predicts that 
as d — > oo the correlation energy vanishes, leaving be- 
hind only the mean-field contribution to the total en- 
ergy Since very little is known about the dimension de- 
pendence of density functionals we cannot say at present 
whether the existence of such a law is a mere coinci- 
dence, a particular property of the Heisenberg model, or 
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a general phenomenon, but we hope that our observation 
of dimensional scaling stimulates further research along 
these lines. 

One practical use that can be made of Eq. (41) is 
to convert an approximate functional obtained for some 
value of d into one for another dimensionality. Counter- 
parts to this property for other Hamiltonians would be 
interesting not only for ab initio calculations (in which 
many results are known for d = 3, but much less in d = 2 
or d = l), 30 but also in the case of the Hubbard model, in 
which the LDA functional is known only for d = 1. 20-22 

IV. CORRELATION ENERGY OF THE 
ANTIFERROMAGNETIC HEISENBERG MODEL 

In this section we apply the results obtained above to 
a study of the correlation energy of the antiferromagnctic 
Heisenberg model. In Sec. IV A we compare the homo- 
geneous expressions of Sec. II with their mean-field ap- 
proximation, to assess the importance and behaviour of 
the correlation energy. In Sec. IV B we study a physically 
interesting inhomogeneity, a spin-density wave, with the 
LDA functional (36). 



(ii) On similar grounds one would expect that corre- 
lations become less important for larger dimensionality. 
This is confirmed both by the main figure and the inset, 
showing that correlations decrease in absolute size and 
relative to the mean-field energy as d increases. The way 
the ^-dependence approaches the classical limit is thus 
qualitatively different from the way the spin-dependence 
does. 

(iii) The mean-field energy is not reliable for any di- 
mensionality d < 3 and S < 5, leading to errors that can 
be larger than J in absolute size and larger than 50% in 
relative terms. This observation puts tight limits on the 
reliability of this rather widely used approximation. 

(iv) The improved treatment of correlations in d = 1, 
represented by the curves labeled 'DMRG fit', does not 
invalidate conclusions (i), (ii), and (iii), obtained on the 
basis of the spin- wave expressions (11) to (13). However, 
both the main figure and the inset of Fig. (4) show that 
it enhances the importance of correlations with respect 
to the mean-field values, as compared to the simpler ex- 
pressions. 

B. Example of an inhomogeneous system: a linear 
spin-density wave 



A. Homogeneous system 

As in the previous section we define the correlation 
energy as the difference between the total ground-state 
energy and its mean-field approximation. For a homo- 
geneous system on a linear, square or cubic lattice the 
latter yields 

E™ F (S,J,d) = -JdNS 2 (42) 

and thus eo(S') = —dS 2 . In the inset of Fig. 4 we plot 
the difference between this value and the total energy 
expressions Eqs. (11), (12) and (13). For one dimension 
we also plot the difference between the mean-field result 
and the more precise expression (15). These differences 
represent the absolute size of the correlation energy. In 
the main part of Fig. 4 we display the relative size of the 
correlation energy as compared to the mean-field energy. 
Several conclusions can be drawn from inspection of these 
curves: 

(i) The inset of Fig 4 shows that the absolute size of the 
correlation energy increases towards larger spins. This 
seems counterintuitive, because larger spins should more 
closely mimick the classical limit, in which there are no 
quantum fluctuations and the mean-field approximation 
becomes exact in the ground state. However, as shown 
in the main figure, the relative weight of the correlation 
energy as compared to the mean-field energy decreases 
towards larger spins. Interestingly, the naive expectation 
that correlations should become less important near the 
classical limit is thus only true in relative terms, but not 
in absolute ones. 



As an example of a truly inhomogeneous situation, to 
which our LDA functional (35) can be applied, we now 
consider a simple but physically interesting inhomogene- 
ity, namely a spin-density wave (SDW) imposed by an 
external field on a chain with antiferromagnctic coupling. 

We model the SDW state by taking 

S„ = Su x cos<p n , (43) 

where <j) n — 2ir(n — 1)/A and u x denotes the unit vector 
in the x direction. This choice describes a linear SDW 
of amplitude S and wave length A, polarized along the 
x-direction. (The lattice is taken to be a chain along the 
^-direction.) The corresponding mean-field energy is 

JV-l N 

£^ F [Si] = J S„ • S n+ i + S„ • B n 

n— 1 n— 1 

= Ey F [Si]+E B [Si], (44) 

where N is the number of lattice sites, and B„ is a 
magnetic field that can be thought of as either exter- 
nally applied [thus forcing the system into a state with 
spin distribution (43)] or generated self-consistently, or 
a combination of both. The LDA approximation for the 
ground-state energy is, on the other hand, 

E& DA [S i ]=E» IF [S i ]+EZ DA [S i ] 
= Ef F [S t }+E B [S l } + J2^ FM (S)\s^ 

71=1 

^ DA [S,]+£ B [S 4 ], (45) 
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where we use, for simplicity, Eq. (36) for E^ DA [Si\. 

For the given spin distribution (43) we now compare 
the predictions of the mean-field and LDA expressions 
for the interaction energy Ej. Note that this is not a 
self-consistent calculation, but a comparison of the two 
expressions (44) and (45) for a fixed distribution speci- 
fied by (43). In Fig. 5 we plot Ef F , Ej DA and E^ DA 
as functions of A, the wave length of the SDW. Since the 
term Eb, arising from the magnetic field, is the same 
in both approximations we only display the interaction 
energy Ej. The presence of an external field with sym- 
metry different from the one of the ground-state of the 
unperturbed antiferromagnetic system gives rise to a rich 
physics. Particularly, we note the following points: 

(i) Addition of the correlation energy to the mean-field 
result lowers the total energy considerably. This is, qual- 
itatively, the same behaviour we found previously in ho- 
mogeneous systems and illustrates again the importance 
of going beyond the mean-field expression for the energy. 

(ii) As A — > oo the SDW approaches a ferromagnetic 
spin configuration. Since the unperturbed homogeneous 
model is antiferromagnetic this state is energetically un- 
favored, and corresponds to a maximum of the E versus 
A curves. In the opposite limit the spatial modulation of 
the SDW can take local advantage of the AFM tendency 
of the underlying homogeneous system. This leads to a 
lowering of the energy as A — > 0. 

(iii) Once A is larger than approximately 10 lattice con- 
stants the E versus A curves saturate. We interpret this 
in terms of the correlation length of the antiferromagnetic 
model by noting that once the SDW modulation takes 
place on a scale larger than a few correlation lengths, the 
system will be relatively insensitive to further approxi- 
mation to the ferromagnetic state. 

(iv) The overall downshift of the LDA curve compared 
to the mean-field one implies that the ferromagnetic con- 
figuration is energetically less unfavorable in the former 
approximation than in the latter. Physically this is rea- 
sonable, because the correlations accounted for by E c 
break up the rigid AFM pattern of the Neel state found 
in the mean-field approximation, and replace it by a com- 
plex ground-state involving spins along all three direc- 
tions in space. 

V. SUMMARY AND OUTLOOK 

Density-functional theory is commonly applied to the 
ab initio Hamiltonian, in the context of electronic struc- 
ture calculations of molecules or solids. 16-18 Applications 
to model Hamiltonians are rare, although they can be 
useful both for the analysis of these models in the pres- 
ence of inhomogeneity (broken translational invariance) , 
and for further development of DFT. The present work 
on the Heisenberg model serves to exemplify these two 
complementary aspects. 

The Heisenberg Hamiltonian is much simpler than 



the Hubbard Hamiltonian, which describes the charge- 
degrees of freedom in addition to the spin ones, or 
than the ab initio Hamiltonian involving the long-range 
Coulomb interaction between charges in a real crystal. 
The Heisenberg model may therefore be considered a 
simplified environment in which concepts and methods of 
DFT can be tested and analysed. An interesting example 
is the study of the interaction and dimension dependence 
of xc functional, about little is known in the ab initio 
case. 

On the other hand, the LDA for inhomogeneous 
Heisenberg models is no more complicated, formally, than 
the mean-field approximation, but it accounts by con- 
struction for essential correlation effects missed by the 
latter. Analysis of the correlation energy of homogeneous 
and inhomogeneous Heisenberg models, as function of 
spin and dimensionality, illustrates that such effects are 
crucial for a quantitative description of the ground-state. 
The above LDA functional may thus be useful in stud- 
ies of the behaviour of the Heisenberg model in spatially 
varying external fields or in the presence of internal mag- 
netic order that breaks translational symmetry. 

The simple LDA-type correlation functional based on 
Eqs. (35) and (11) is seen to produce significant quantita- 
tive change and qualitative improvement over the mean- 
field approximation, at very little extra computational 
cost. This observation encourages us to envisage more 
complex applications of this functional, such as to im- 
purity states in the Heisenberg model. An extension of 
the present work to a study of the thermodynamics of 
inhomogeneous Heisenberg models (employing the T > 
formulation of DFT 31 ) is also planned for the future. 
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FIG. 1. Per-site ground-state energy of the homogeneous 
ferromagnetic Heisenberg model in one, two and three dimen- 



FIG. 2. Per-site ground-state energy of the homogeneous 
antiferromagnetic Heisenberg model in one, two and three di- 
mensions. Dashed lines: interval according to Eq. (8). Full 
lines: spin-wave theory expressions (11), (12) and (13), re- 
spectively. Full circles (in one dimension only): DMRG data 
of Ref. 13. The inset is a zoom into the low spin region of the 
main figure. 



FIG. 3. Same as in Fig. 2, but on a reduced scale. On the 
vertical axis the expression (11) has been subtracted to dis- 
play only the corrections to that result. Full circles: DMRG 
data of Ref. 13, dashed line: fit proposed in that reference, 
present Eq. (14), full line: fit proposed here, Eq. (15). 



FIG. 4. Correlation energy of the antiferromagnetic homo- 
geneous Heisenberg model in one, two and three dimensions, 
as a function of spin. The curves labeled d — 1,2,3 are ob- 
tained by subtracting the mean-field energy from the expres- 
sions (11) - (13), and the one labeled l d = 1 DMRG- fit' is 
obtained by using our fit Eq. (15) instead of Eq. (11). Main 
figure: correlation energy relative to mean-field energy. Inset: 
absolute size of correlation energy. 



FIG. 5. Interaction energy Ej of a spin-density wave im- 
posed on an antiferromagnetic chain, as a function of wave 
length of the spin modulation. Upper curve: mean-field ap- 
proximation, defined by Eq. (44). Lower curve: local-density 
approximation for the correlation energy, as given by Eq. (36). 
Middle curve: LDA result for interaction energy, obtained as 
sum of the two other curves. 
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